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Abstract 

In this paper, we present new results relating the numerical range of a matrix A with 
generalized Levinger transformation C(A, a, 0) = chHa + (3Sa, where Ha and Sa, are, 
respectively the Hermitian and skew-hermitian parts of A. Using these results, we, then 
derive expressions for eigenvalues and eigenvectors of the perturbed matrix A + C(E, a, 0), 
for a fixed matrix E and a, (3 are real parameters. 

Keywords : Numerical range, Generalized inverses, Perturbation theory, Eigenvalues-eigenvectors, 
Control and Systems Theory, Sensitivity. 

AMS Subject Classifications : 15A60, 15A09, 47A55, 65F15, 93B60, 93B35. 

1 Introduction 

Let A4 n (C) (Ai n (M)) be the algebra of all n x n complex (real) matrices, and let A G 
A4 n (C). The numerical range of A, also known as the field of values jB], is the set 

NR[A] = {x*Ax G C : x G C n with x*x = 1}. 

The numerical range NR [A] is a compact and convex subset of C, that contains the spectrum 

cr(A) of A. If A G o~(A) PidNR[A] with multiplicity m, then A is a normal eigenvalue, i.e., A 

is unitarily similar to XI m (BB, with A ^ o~(B). Clearly, when A is normal, then its eigenvalues 

are normal and Ai?[j4] = Co{a(A)}, where Co{ • } denotes the convex hull of the set. For 

any A G A4 n (C), if we write A = Ha + Sa, where 

A + A* A -A* 
Ha = 7, and da = , 
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are the Hermitian and the skew-hermitian parts of A respectively, then 

Re NR[A] = NR[H A ] and i Im NR[A] = NR[ S A ] ■ 

Moreover, if P is any n x m matrix with n>m and P*P = I m , then 

NR[P*AP] C NR[A], 

and the equality holds only for m = n. 

Given a matrix A G A4 n (C), we define "generalized Levinger transformation " of A as the 
double parametrized family of matrices 

C(A,a,P) = aH A + (3S A = ^^-A + ^^-A*, with a,f3 G R. (1) 
In (dJ, for a = 1 and /3 = 2t - 1, t 6 R, we have 

£(A, 1, 2i - 1) = tA + (1 - t)A*, 
i.e., jC(A, l,2t — 1) is just the ordinary Levinger's transformation, which has been studied in 

nana nam]. 

Moreover, for a = 2t — 1, t £ R and (3 = 1, we have 

£(A,2f-l,l) = tA + (t-l)A*. (2) 

The equation is a different formulation of Levinger transformation, where in C(A, 2t — 1,1) 
the difference of coefficients t and t — 1 of matrices is equal to unity (skew convex expression) . 
Note, that 

£(L4, 1, 2t - 1) = iC(A, 2t - 1, 1). 
Clearly, for every a,j3 6 I, we have from (JTJ) 

Hc{A,a,p) = CiH A , Sc{A,a,p) = &S A - 

Hence, 

ReiVi?[£(A,a,/5)] = NR[H c{Aa ^] = aNR[H A ] = aReNR[A] 
ilmNR[C{A,a,(3)} = NR[S C (A,a,p)] = ifiTm NR[A], 

and consequently 

NR[£(A,a,/3)] = {ax + if3y : x,y G R , with x + iy G iVi?[A] }. (3) 
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Moreover, the boundary of NR[£(A,a, (3)] is given by 



dNR[£(A,a, 



{ax + i(3y : x,y £ R, with x + iy G dNR[A] }. 



(4) 



For A G M n (R), since NR[A] = NR{A T ], NR[C(A,a, (3)] is symmetric with respect to 
the real axis and 

NR[£(A,a,P)] = NR[£{A*,a,-f3)], 

the domain of (3 can be reduced to [0, +oo). Additionally, if < (3± < /?2, then due to ([3]), 
@ and the symmetry of the numerical ranges with respect to the real axis, we have (in some 
sense) a vertical dilation, i.e., for z\ € dNR[£(A,a, fii)], Z2 £ dN R[£(A, a, fa)], holds 



•j - \l a 2 {x*H A x) 2 + @1{x*Sax) 2 < \l a 2 {x*H A x) 2 + (3%(x*S A x) 2 = \z 2 \ 
and consequently 

NR[£(A,a,@i)} C ATi2[£(A,a,#j)]. 



Example. For ^4 



the numerical ranges of A and a, /3) are illustrated 



13 4 
2 7-6 
-13 5 

in the following figures. On the right, due to ([3|), the vertical dilation is presented only when 
a is fixed (here a = 0.4) and j3 is altered (0.5 < (3 < 1.2), reminding that, this property holds 
only for the ordinary Levinger's transformation. Otherwise, NR[£(A, a, (3)] is moved as it is 
shown on the left figure, for the values a = —0.9, (3 = 0.8, a = 1.3, [3 = 0.6 and a = 1.4, (3 = 1.3. 
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The paper is devoted to the study of the generalized Levinger transformation of a matrix. 
Specifically, we establish an interesting relationship between the numerical range of a matrix A 
and its generalized Levinger transformation. This relationship is then used to obtain results on 
the eigenvalue and eigenvector of the perturbed matrix of the form A + C(E, a, (3), where E is 
fixed and a, and (3 are small parameters. 

Our motivation for such study comes from the fact that a great deal of effort has been made 
in the literature to establish bounds on the eigenvalues of a perturbed matrix. For results on 
this topic, see [13] and the well-known books on linear and numerical linear algebra by Datta 
[3], Stewart and Sun [T3], Kato [7J, Lancaster and Tismenetsky [8], and Bhatia [2]. 

This paper is divided in two parts. The first part contains geometric properties of numer- 
ical range of C(A, a, (3). Also bounds are given for real and imaginary parts of eigenvalues of 
C{A, a, (3). This provides us a framework to study variation of the spectrum of C(A,a, (3). In 
the second part, we use the Levinger transformation for a fixed matrix E as a perturbation 
matrix, whose activity on a matrix A depends only on the parameters a and (3. First, we 
formulate a necessary and sufficient condition for a normal matrix to remain normal, under 
a perturbation by a symmetric and rank one matrix. Next, we present an approximation of 
a perturbed eigenpair of a diagonalizable matrix A using two parameters, which generalizes a 
known result in [13], p. 183], where the eigenvector of perturbed eigenvalue is not mentioned 
and even the perturbed eigenvector is investigated in [4, p. 431], without giving further details 
for the perturbation of the corresponding eigenvalue. Further we simplify these formulae using 
the notion of generalized inverse extending the corresponding result in [S]. As an application, 
we present a sufficient condition such that the perturbed eigenpair is first order approximation 
of the corresponding simple eigenpair of initial matrix A, and we give two numerical examples 
to illustrate our results. 

2 Geometric properties 

Proposition 1 a. Let A £ M n (C). The image of a line segment €a € ./V.RL4] by the Levinger 
transformation, is a line segment ec £ N R[C{A,a, (3)\. 

h. If A £ M n (M) and NR[A] is an ellipse, then NR[C(A, a, (3)], for is also an ellipse. 

Proof, a. For any x\ + iyx,x 2 + iy 2 £ ./V\RL4], and t G [0, 1], we observe that 

a[(l - t)x\ + tx 2 ] + /3[(1 - t)iyi + tiy 2 ] = (1 - t)(ax\ + i/3yi) + t(ax 2 + i/?y 2 ), 

where ax\ + \(3y\ and ax 2 + i(3y 2 he in the convex set NR[C(A,a,f3)]. Hence, the proof of a 
follows readily. 
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b. Consider that NR[A] = {x + iy : ^ + < 1, with c > k > }. For a / 0, and 
(3 > 0, changing the variables x = a X, y = (3~ 1 Y, by ([3]) we have that the boundary 
dNR[C(A, a, /?)] is the ellipse + = 1. The foci are on the real axis, when a £ 

M\[ — -@— , -^r-], otherwise, they lie on the imaginary axis. □ 

Proposition 2 Let A £ M n (C), B £ M m (C), and let NR[B] be a polygon circumscribed to 
NR[A\. Then the geometric relationship of numerical ranges NR[C(A, a, (3)] and NR[C(B, a, (3)} 
remains the same, for a,/3eR - {0}. 

Proof. By Proposition QJi and (j3|), NR[C(B, a, (3)] is a convex polygon. Moreover, it is 
easy to see that there do not exist other common boundary points of NR[C(A, a, (3)] and 
NR[C(B, a, /?)], except those, which correspond to z £ dNR[A] n dNR[B]. □ 

Proposition 3 Let A £ M n (C), and a, (3 £ K - {0}. 

a. If A is a normal matrix, then dNR[C(A, a, (3)} is a k— polygon, as dNR[A\. 

b. NR[C(A,a,P)] nR = Kn NR[aA]. 

Proof, a. The proof of this part follows readily by Proposition [Th and from the observation 
that A is unitarily similar to diag{ x\ + iyi, X2 + iy2, . . . ,x n + iy n } if and only if C(A, a, (3) is 
unitarily similar to diag{ ax\ + i/3yi, ax2 + i(3y2, • • • , ctx n + i(3yn}- 

b. Clearly, if ax* Ax £ NR[aA] n E, then x*£(A, a, (3)x = x* Ax + x*A*x = ax* Ax, 
concluding that NR[aA] flic NR[C(A, a, (3)] nt. Moreover, for M\{0}, because 

o^4 = a, (3) - C*(A, a, (3), 

in an analogous way, we obtain NR[C(A, a, (3)] D R C iV\R[a<yl] fll. □ 



Remark 

The eigenvalue A S cr(^4) is normal if and only if A£ is a normal eigenvalue of C(A, a, (3). 
In fact, from and the relationship U*AU = \I m © i?, where A ^ &{B), it follows that 
a, /?)[/ = X c I m © £(5, a, /?). 

In the following proposition we present a compression of NR[C(A,a, (3)], when A is a 
normal matrix, based on our results in pQ. 

Proposition 4 Lei A £ M n (C) 6e a normal matrix and let the polygon < Af^j , A(^) 2 , . . . , A(^) fc > 
be the numerical range of A. If xj is a corresponding eigenvector of A^)., j = 1,2, ...,k, 
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and v = Y2j=i v j x j * s a un ^ vector, denoting by E = span{ v } and Eyy the orthogonal 
complement of E with respect of W = span{ x\ , X2 , ■ ■ ■ , Xk }, then 

NR[P*£(A, a, 0)P] = NR[C(P*AP, a, /?)] C < A (£)l , A (£)a , . . . , \ c)k >, 



where P 



[w\ W2 ■■■ iOjfe_i], and w\ , W2 , • • • , "Wfc-i is an orthonormal basis of E^. 

(£)i , A(£) 2 , . . . , A (£)fe 



Moreover, dNR[P*C(A,a, /3)P] is tangential to the edges of the polygon < X(c)i > A 
at the points 



> 



M(£)t = aRe^ {A)T + if3Imii {A)T , (r = 1, ...,/V) 



(5) 



where 



M(A)t 



br+il 2 A(^) T + |-u r | 2 A (A ) t+1 



|iV+i| 2 + |u T | 2 



(r = !,...,&- 1 



M 2 A (A)fc + K| 2 A (A)l 



Proof. By the equation (JTJ) we have C(P*AP, a, (3) = P*C(A, a, /3)P, and it is evident the 
equality of the numerical ranges. Moreover, C(A, a, (3) is normal, and it is known in [I] that, 
dN R[P* C(A, a, P)P] tangents to the edges of the polygon < X(ch , \c) 2 >■ ■ ■ j \jC) k > a t the 
points 

K +1 | 2 A (£)r + K| 2 A (£)t+1 |ui| 2 A (£)fe +\v k \ 2 A (£)l 
r- i2 T i T , i2 (r - 1) , M(£)fc - — 

I <Jt+1 | t | u T | 



kll 2 + kfcl 2 



Since, the eigenvalues of A and C(A,a,f3) are related by \c)r = a R e A(A) T + fc/3ImA( A \ T , the 
equation © is verified. □ 



In the following applying the results of Rojo and Soto in [12] to C(A,a,(3), one obtains 
bounds for the real and the imaginary parts of the eigenvalues of C(A, a,/3). 

Theorem 1 Let the matrix A 6 M n (M) and Xj 6 a(C(A, a, /?)). Then for each Xj we have 

1/2 



ReXj 



and 



atr{H A ) 



n 



< a 



\ImXj\ < \/3\ 



n — 1 



71—1 



72. 



I-zTaIIf 



ll^lll 



3(o?\\H a \\ 2 f + IP\\Sa\\ 2 f ) 



22 



a 2 ||5 A flA-fl- A 5 A 



|2 



3(a 2 ||^ A || 2 , + /? 2 || > S A || 2 7 ) 



1/2 



(6) 



(7) 



where 



denotes the Frobenius norm. 
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Proof. Observe that tr(£(A,a, P)) = tr(aH A + PS A ) = atr(H A ), and 

\\£(A, a, P)£ T (A, a, (3) - £ T (A, a, P)£(A, a, /3)\\ F 

= \\(aH A + (3S A )(aH T A + /3S%) - {aH T A + pS T A )(aH A + (5S A )\\ F 

= 2\af3\\\SAH A -H A S A \\ F . (8) 

Since 

tr{H A S T A ) = —tr(H A S A ) = -tr ^ + A?)(A - A^ = _X_ ^ _ ^ [( ^ T)2]) = Q 
we have 

\\£(A,a,P)f F = tr[£(A,a,/3)£ T (A,a,f3)]=tr[(aH A + f3S A )(aHl + pSl)] 

= a 2 \\H A \\ 2 F + (3 2 \\S A \\ F (9) 

and 

tr[C 2 {A,a,P)\ =tr[{aH A + pS A ){aH A + pS A )] = a 2 \\H A \\ 2 F - p 2 \\S A \\ 2 F , (10) 
and consequently by Q and (fTUj) 

\\C(A, a, P)\\ 2 F + tr(C 2 (A, a, 0)) = 2a 2 \\H A \\ 2 Fl 

(11) 

\\£(A,a,p)f F - tr(£ 2 (A,a,P)) = 2[3 2 \\S A \\ F . 
Therefore, if we substitute ([8]), ©, and (jlip in the relationships of Theorem 7 in |12j . we have 
tr(£(A,a,/J)) 



ReAj 



and 



< 



' n-1 / \\£(A,a,P)\\ F + tr(£ 2 (A,a,P)) _ v(£(A,a,P))* _ [tr(£(A,a, P))]< 
n V 2 12||£(A,a,/J)||J. 



|ImA,| < J«Z1 (\\£(A,a,PW F ~ tr(£ 2 (A,a,P)) - 



where v(£(A,a,p)) = \\£(A,a, P)£ T (A,a, P) - £ T (A, a, p)£(A,a, P)\\ F , thus we obtain the 
bounds for the real and the imaginary part for each eigenvalue of £(A, a, P) in © and ©. □ 
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3 Application to perturbation theory 

The question " how close is a matrix M to being normal " , it is known that it is evaluated by 
the normality distance \\AA T — A T A\\ p of p norm. Since various matrix norms are equivalent, 
using the Frobenius norm we have: 

Proposition 5 Let N £ Ai n (M.) be a normal matrix and for a nonzero vector x £ M. n , 
let E = xx T . If x is not eigenvector of N corresponding to a real eigenvalue, the matrix 
M = N + E is normal if and only if N is symmetric. 

Proof. For the symmetric matrix E = xx T , clearly E 2 = \\x\\ 2 E, and the normality distance 
of M is equal to 

\\MM T - M T M\\ F = \\(N + E)(N T + E) - (N T + E){N + E)\\ F 

= \\NE + EN T - N T E - EN\\ F = \\R + R T \\ F , 

where R = NE - EN. Since 

tr(R 2 ) = 2 [tr(NENE) - tr(EN 2 E)] = 2 [{x T Nxf - ||x|| 2 (x T iV 2 x)] 
tr{RR T ) = 2 [ \\x\\ 2 tr(NEN T ) - tr(ENEN T )] = 2 [||x|| 2 (x T iV T iVx) - (x T Nx) 2 } 

and tr [(R T ) 2 ] = tr(R 2 ), we have : 

\\R + i? T || 2 = tr [(R + R T ) 2 ] = tr(R 2 ) + tr [(R T ) 2 } + 2tr{RR T ) = A\\x\\ 2 x T (N T N - N 2 )x. 

Hence, the matrix M is normal if and only if N T N = N 2 . This equation is equivalent to 

DD = L> 2 , (12) 



where D is diagonal and unitary similar to N, i.e., N = UDU*. Thus, by (|12j) . D is real and 
N is symmetric, since it is unitary similar to a real diagonal matrix. □ 

Corollary 1 Let N = diag(Ni, N2, • • • , N T ) G M n (R) be a normal matrix and for a nonzero 
vector x = xT . . . £ M n , with all Xj 7^ 0, let E = xx T . If Xj is not eigenvector of 

Nj, (j = 1, ... ,t), corresponding to a real eigenvalue, the matrix M = N + E is normal if and 
only if N is symmetric. 

It is worth notice that, the result in Proposition [5] is combined by the special form of E, 
but it is interesting to look at more general perturbations, investigating how main properties of 
./V are influenced. For this, we consider E £ A4 n (R), since 

\\MM T - M T M\\ F = \\2H [N ^ T] + EE T - E T E\\ F , 
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where [N, E T ] = NE T — E T N, we conclude that, the normality distance of M is related to 
the normality distance of E. Hence, an outlet is to investigate if some properties of perturbed 
normal matrices remain. 

Let the matrix A £ Ai n (C) be diagonalizable (keeping the property of N) and E 6 Ai n (M) 
be fixed, without giving any attention to \\E\\. Consider the matrix 



where a, (3 G M are small enough varying parameters. Clearly in (j 1 3 j> . M a g is continuous 
differentiable and the Hermitian and the skew-hermitian parts of E influence independently the 
matrix A. Especially, when A is normal, He and Se alter Ha and Sa separately. 

Denote by X a ^ an eigenvalue of M aj p in (fT3j) and by v a ^ and u a> /3, the corresponding right 
and left eigenvectors, i.e., (M a ^ — X a ^I)v a ^ = 0, p(M a ^ — X a ^I) = 0. Since the coefficients 
of characteristic polynomial det(XI — M a p) are polynomials of two variables a, (3 and A Qj/ g is 
continuous function of these coefficients, for a = (3 = the perturbed eigenvalue X a ^ is equal 
to a semisimple eigenvalue Aj of A, and the eigenvectors are : v a ,p = Vi, = cui, where Vi 
and uji are the right and left eigenvectors of Aj for the matrix A. We remind the readers that 
an eigenvalue is called semisimple, when it is a simple root of the minimal polynomial of matrix. 
Moreover, A Q)( g and v a ^, u a ,p are continuous functions of a,/3 and partial differentiable, but 
might have rather singularities on total differentiability [7, p. 116]. For further details we refer to 
[7] and [SJ Ch. 11]. We will now give a result on the sensitivity of eigenvalues and eigenvectors 
of perturbed matrix M a ^ in (|13p in the neighborhood of A, in relation with the remaining 
eigenvalues and eigenvectors. 

Theorem 2 Let the matrix A S M n (C) be diagonalizable and let Vj and ujj be the right and 
left eigenvectors of A corresponding to Xj £ cr(A). If the eigenpair (A a g , v a s) has continuous 
second order partial derivatives in the neighborhood of Aj and Vi, then : 



M a3 = A + C(E, a,(3)=A + aH E + 0S E , 



(13) 




Xi + 



uj*£(E, a, P)vi 



(io*C(E,a,f3)v k )(uj* k C(E, a ,p)v i ) 
(Xi — Afc) SjSfe 



(14) 



k^i 



v k ujlC(E,a,P)vi 





(u* C(E, a, g) y k )(u* £{E, a, (3) Vj) 
(Xi — Aft)(Aj — Xj) s k Sj 



-E 



(lo* £{E, a, /3) v l )(u* £{E, a, g) Vi ) 
(Xi — Aj) 2 SiSj 



v j + 0(a 3 ,p 3 ) 



(15) 



where S£ = u>gV£. 
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Proof. The partial derivatives of the equation (M a ^ — \ a ,pl) Va,/3 = 0, with respect to a, (3, 
are 



(16) 

8(3 ' J ' ^" ,,J d/T = ° 

Multiplying these by cj*^, since u* a ^M a ^ = \ a ^io* a ^, we have 

</3 (he ~ A v atf ) = 0, u* j/3 (s E ~ ^jf 1 ^ v <*,0 = 0. 

For (a, /3) — ► (0, 0) the expressions given above 

d \*,P)=o = lim d\ a:f3 = uj*H E Vj d\ a ,i3)=o = nm d*a£ = U*S E Vi 
da (Qt,/3)-»(o,o) 5a w*^ ' 9/3 (a,/3)->-(o,o) 9/3 w*t>i 

and then, the first differential dX a ^ is equal to 

„ d \a,/3)=o , d\a,/3)=o uj*H E Vi uj*S E Vi Lu* £(E , a, /3)vi 
dX a ,f3 = a 5 h p ^ = a ; 1- P — — 



da d(3 uj*Vi oj*Vi uj*Vi 

Moreover, the first equality in (fT6j) for (a, (3) — > (0, 0) gives 

Since ^4 is diagonalizable, we can write = 5^J* =1 c^i^, and so the last equality can be 

written as 



Furthermore, multiplying the above equality by the left eigenvector uJk of A, and using the 
orthogonality of u>k and Vi (k ^ i), we have 

c fc = u k H E v % — for fc ^ 



and consequently, 



&W)=o _ V- w* k H E Vi , . 
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Similarly, by the second equality in (fT6|) . we obtain (Aj — X k )oj* k — -^jw ° = ^lS E Vi, and thus 

frW)=o _y^ u* k S E Vi , . 

9/3 ^(A.-A,)^^ ^ 

Hence, the differential cfo aj( g can be computed as 

LO* k H E Vi sr^ u* k S E Vi _sr^ v k u* k C(E 1 a,f3)v i 



^(Ai-Afc)^ ^ (Ai-A fe )o;*y fc ^ (Ai - A fc ) wjE 



Now, the partial derivatives of the equations in (|16p with respect to a, /?, are 

( d\ a 3 \ 8v a 3 , 8 2 v a 3 8 2 X a 3 

2 { He --8^')^ L + < M «." - ^ ' > "S? 2 - "ft? 2 = 



2 ( - — I j — + ( M a , 3 - X a ,3l )-^—^- v a ,3 = (20) 



dad/3 \ da J 8(3 \ 8(3 J da da 8(3 

Multiplying these expressions by ^ and substituting 9v{o q^ =0 , dV{a ^=° from (fT8|) and JTH 
for (a, (3) — > (0,0), and noting that = 0, we obtain 

<9 2 A (Qj/3)=0 2 / t ^W)=o\ *(a,)3)=o 2 ^{ujlH EVl ){uj*H E v k 



9a 2 V 1 ~' ^« 7 ^« jbi (Ai-A fc )u;|f fe 

<9 2 A(a,/?)=o 2 / ^A (ai/3)=0 ^ dv( ai p )=Q 2 x - (^S E Vi)(u}iS E Vk) , on 



d/3 2 wj-t;. V 1 9/3 y 5/3 (A.-Afc)^ 

9 2 A( Q ,/3)=o 1 f *„ <^K/3)=o 9u( a)i 9)=o 

ftt; h Wjiis 



<9ad/3 V d/3 1 da 

_ 1 ^fc )( U* k S E Vj) + { U* k H E Vj )(u*S E v k 

uj*vi ^ (Aj - X k ) uj* k v k 

Therefore, the second differential d 2 A Qj/ g is equal to 
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2a 2 (LO* k H E Vi)(uj*H E v k ) 2a (3 {u*H E v k )(uJ* k SE Vj) + ( i^flg ^ ) {uj*S e Vk ) 
u * v ii^i (Ai-A fc )u;*^ u*Vi j£ (Xi- X k )uj* k v k 

2(3 2 ^ (u* k S E Vi){u;tS E v k ) 



(Xi - X k ) uj* k v k 

2a_ (^^^(^(^.a,/?)^) 2ff ^{u*S E v k )(w* a, /?) ) 



and by 



U i V i ^ ( A * - A *0 ^fcf fc ^ (Aj - X k ) C^fc 

2 K «, g) ^ ) K a, g) v t ) 
(Ai-A fe )a;*^ 



A a ,/3 = Ai + dA Qj/3 + yd 2 A a>/3 + (a 3 , /3 3 ) 



we receive (|14p . whereas we have declared sg = lo»V£. 

Multiplying the first of ([20]) by coj, due to ufa = (j / i), for (a,/3) — > (0,0), we obtain 



5a 

Substituting the formulae of d \%> = ° , ^f^ from (HTJ) and CED, since cu*u fc = (j ^ k), 
we take 



» dy,fl=o _ / V HH EVl )(u;*H E v k ) \ 2{u* i H E v i )(u* j H E v i ) 
U}j da* 1 ^ (A, - A fc )(A< - A,) u%v k J u*v t (X t - A,) 2 ' 3 * % 



and then 

52 



g \*,f})=o _ I V {^lH E Vi)(u:*H E v k ) \ (u,;H E Vi)(u]H E Vi) 

a* 2 h\h - A *)( A » - A i) K^) ) Wj ^ - A ^ 2 K t*)(«;«i) Url j 



Similarly, the last two expressions of (|20p lead to 



g«M=0 = y ( {.^EViX^jSEVk) \ (uj*S E Vj)(uJ*S E Vi) 

= f£ [t^i ^ ~ X ^ ~ A i)K^) H^-) J Ui ^ ( A * - A . ) 2 « u i) Vj 

(23) 

&V{a,P)=0 V" | V ^fc 5 ^ v i)(Vj H E Vk) + K^fi Vj)(UjS E Vk) \ 

- jg I 2^ (Aj - A*)(A, - A,-)K^) (a,**,,-) J Uj 



12 



\{uj*SEVi){iJ*H E Vi) + (u*H E Vi){Ld*S E Vi) 



Thus, by ([22]) and ([23]) we take 

\- \- {uj* H E v k ){uj* k C{E,a,P)v i ) (uj* S E v k )(iv* k £(E,a,p)vi) 

2a v (w* H E Vj)(u>* C(E, a, (3) Vj) _ 2$_ (tj S E Vj)(uj* £(E,a, f3) Vj) 



(A i -A fe )(A l -A J )K^)(^) Uj 

u,*^ (Aj — Ay) 2 (a>|t!j) 



and by 



we obtain the claimed equality ([15]) . □ 



V = Vi + dv a ^ + ^-d 2 v a ^ + (a 3 , /3 3 ) 



The simplified presentation of partial differential formulae (]17p and (|21[) of A Qj/ g and (118p . 
(HH), ([22]) and ([23]) of u Q / g for a = /? = are independent results on those, which were obtained 
earlier by Chu in [3]. Chu has follow different methodology considering that Aj is simple 
eigenvalue and an additional normalized condition that uj*V{ = 1, and even Chu's formulations 
of the partial derivatives depend on the invertibility of a matrix and the eigenvectors t>j, Up 
Furthermore, no results on the perturbation of the eigenpairs X a ,p, and v a s are given in [3]. 

In the following we present a lemma, which will contribute in the approximation formulae 
CE1 and ([S]). 

Lemma 1 Let the matrix A E M n (C) be diagonalizable and YJ, Wi be matrices whose 
columns t>j and rows to* respectively are the corresponding right and left eigenvectors of A, for 
Xi G cr(A). A generalized inverse of (A — Aj/) M , /i G N, is defined by 

[(a - xar] + = y, 7T Vk f„ ; *k = ufa. (24) 

(Xk-XiYsk 
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Proof. It is evident that (A — ^'-a. = v k oj^., and then 

V ; (A fe - XiY v ; 1 ; (A fc - \)" 

= (A -Ail)" 77 r^— T = ••• = v k u k . 

Since Yl = ^ = AjW* and for fc ^ i, u* k Vi = 0, we have: 



A' 



and 



(a - xar - YiWi(A - xar = (a- xar - (w t A - wo (a - w- 1 = (a - a< 



V M _ A /^ V 



V \(I - Y W-) = V 

^(x k -x t y Sk ) { 1 ij Z. (Afc -A0^ fc - 

□ 

In Lemma [TJ if A is normal, then v k = u) k , and [(A — Aj/) M ] + is Hermitian. In this case, 
we confirm that [(A — Aj/)^] + in (|24p is the Moore-Penrose inverse of (A — Aj/)^. 



Combining Equations (|24p , (|14p , and (|15p , Theorem [2] leads to a generalization of a corre- 
sponding result for simple eigenvalues of a Hermitian matrix, which was presented by M. Fiedler 
in 0. 

Theorem 3 Let the matrix A E A^n(C) fre diagonalizable and Aj 6e a semisimple eigenvalue 
of A with Vi, uji corresponding right and left eigenvectors. If the assumptions for the equations 
(TT^P and U5\) hold, then the following expressions for A aj/ g and v a $ hold: 

X a ,f3 = X i + —i0*C(E,a,p)v i -—L0*C(E, a ,P)(A-X i I) + C(E,a : p)v i + O(a 3 ,f3 3 ), (25) 

&i Si 

v a ,fs = v i -(A-X i I) + C(E 1 a,f3)v i +[(A-X i I) + C(E,a,(3)} 2 Vi 

- — [(A - A,/) 2 ] + C(E, a, (3) v^* C(E, a, P)v l + (a 3 , /3 3 ) . (26) 



14 



Proof. From (|14p and (|24p with [1 = 1, we immediately have 

\ a ,p = A,, + — a, - — uj*£{E, a,p) I V j^^k ) «, £0^+0 (« 3 , /? 3 ) , 

proving (|2"5|) . Also, from ([To]) and (|2~4"|) with /i = 1,2, we have 



E 



' (Afc — Aj) s/c 



-C(E,a,p) n VkU) ,\ £{E,a,/3)vi 



(Aj — Aj) Sj 



(Afc — Aj) Sk 



E 



£(E,a,P)^^£(E,a,P)v l + 0(a 3 ,p 3 ) 



V, 



(A - Xil)+ C(E, a, 0) Vi + {A- C(E, a, (3) £ 



v k uj k 



k^i 



(Afc — Aj) Sfc 



£(E,a,p) Vi 



E 



VjU- 



(Xj - Xi) 2 Sj 



C(E, a, (3) ^ £(E, a, 0)v l + O (a 3 , /3 3 ) 
s.- 



= fj — (A — Aii)+ a, (3) Vi + [(A - Xil) + £(E, a, (3)] v t 

- [{A - XJ) 2 ] + £{E, a, 0) ^ £{E, a, f3) Vi + (a 3 , /3 3 ) , 

Si 

proving (|26p . □ 
In (|25|) and (|26p . if we consider the first-order approximation, then we simply have 



1 



\*,p = \ + —u*£(E, a, (3)vi + (a 2 , /3 2 ) 

Si 



and 



(27) 



v a , p = Vi -(A- XiI) + £(E, a, P) Vi + {a 2 ,p 2 ) 
for a simple eigenvalue Aj. In these cases, 



M a ^v a ^ - X a ^v a ^ = -(A - XJ)(A - XiI) + £(E, a, f3) v { + £(E, a, (3)vi [uj* £{E, a, (3) Vi ] v { 

Si 



+ 



(cj*£(E,a,(3) Vi)I- £{E,a,(3) 



(A-XiI) + £(E,a,(3) Vi 



(A-X l I)(A-X l I) + + V ^ 



s. 



£(E,a,(3) Vi+£(E,a,f3)vi 



+ 



-(oj* £(E, a, p) Vi )I - £{E, a, p) 

Si 



{A - XiI)+£(E, a, p) Vi + O {a 2 ,p 2 ) . 
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Since, 



(A - \iI)(A - XJ) + + ^ = (A- XJ) V 



'. (Afc — Aj) Sfc 



+ 



we have 



-u* C(E, a, (3) va - C{E, a, (3) 



(A - \T)+C{E,a,P)vi + O (a 2 ,/? 2 ) .(28) 



Proposition 6 Let Aj be a simple eigenvalue of diagonalizable matrix A £ M n (C) with right 
and left eigenvectors Vi and u>i. If there exist a, (3 such that C(E,a, (3) Vi £ ker(A — Aj/) + , 
then X a ^ and v a fi in \21^ is an approximation of an eigenpair of M a ^ = A + C(E,a, f3). 

Corollary 2 Let Aj be a simple eigenvalue of normal matrix A £ A4 n (C) with eigenvector V{. 
If there exist a, {3 such that £(E, a, (3) Vi = Vi, then and v a ^ in [27\ ) is an approximation 
of an eigenpair of M a ^ = A + C(E, a, (3). 

Proof. It is well-known that 

ker(A - Xil) + = ker(A - XJ) T . 
Also, since A is normal, uJi = V{. Thus by (I28p . it is implied that 

(A - XJ) + C(E, a, (3) vi = (A - XJ) + v\ = [v* (A - XJ)f = 0. 



□ 



In this case, (|27p is simplified to 



A Q ,/3 = Ai + ^ + O (a 2 , (3 2 ) , v ajP = Vi + G (a 2 , (3 2 ) 

Si 



(29) 



or, X a ,p = (Xi + l) + 0(a 2 ,(3 2 ) 



Vi 



+ O (a 2 ,/? 2 ) for real symmetric matrix A. 



Example 1 Let A 



1 1 
2 
2 

right and left eigenvectors are given by : 



. Then a (A) = { Ai = 1, A2,3 = 2}, and the corresponding 



Vl v 2 v 3 



1 1 
1 
1 



UJl UJ 2 U>3 



1 
1 
-1 1 
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Let the perturbation matrix E 



1 





1 







1 





. We then obtain 


-1 











£(E,a,f3) 



a 
a 
-P 



and M, 



a,/3 



1 + a l+P 
2 + a 
-P 2 



„YJLf ^ /\ 3+ a -V(a-l) 2 -4/3(/3+l) 3+a+V(a-l) 2 -4/?Qg+l) 

Inen, <T(M aij aj = {Ai iaij g = Y ^ > A 2,a,/3 = - § 

^3,a,/3 = 2 + a}, and the corresponding right eigenvectors are: 



l + P Ai a R -I- a 



l + P A 2 , Qi/ 3 - I- a 



1 



Clearly, the eigenvalues and the corresponding eigenvectors of M a ^ are real functions of two 
variables, with continuous partial derivatives for all permissible values of a, p. 
For a = 0.1, P = 0.01, we have a{M a ^) = { 1.1114, 1.9886, 2.1 }, and the corresponding unit 
eigenvectors are given by 



0.9999 0.0113 
Moreover by ([24"D . 

(A - /)+ = [(A - I)Y 



0.7508 0.6606 



1 



1 
10 
1 



Hence, a first-order approximations of the eigenvalue and unit eigenvector of M a p by (f2"T|) are: 

1 T 



Ai = 1 + 0.11 = 1.11, v{ 



1 0.0099 



Also, by ([25]) . ([26]) . the corresponding second-order approximations are equal to 

l T 



By HMD, (A -21)- 



Ai = 1.1112, 

-1 1 





0.9999 0.0111 



, P-2/) 2 ]- 



1 -1 





and substituting these 



expressions in (j25[) and (j26H , using the eigenvectors u 2 , W2 , the second-order approximations of 
the eigenpair of M a ^ are 



A 2 



^2 



0.7505 0.6609 
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Similarly, using the eigenvectors ^3,0^3 , the second-order approximations of the third eigenpair 

of M Qj/ 3 are 



A3 = 2.1, vi 



1 



Notice that, we obtain all eigenpair of M a ^ with preciseness 10 



-3 



Example 2 Let the matrix A 



2 


-2 





-2 


6 


-1 





-1 


2 



and 



V\ v 2 ^3 



1 2 2 
1 -5 
-2 1 1 



with cr(A) = { Ai = 2, A 2 = 1, A 3 = 7}, 



be the corresponding right and left eigenvectors. 



For a = 0.04, (3 = 0.08, and non symmetric £7 



£(£7,0.04,0.08) 



0.84 0.16 -0.08 

0.04 
-0.4 0.8 



and M, 



21 3 -4 
110, then 
-8 20 

2.84 -1.84 -0.08 

-2 6.04 -1 

-0.4 -1 2.8 



a,f3 



We have a(M a> p) = { 3, 1.5886, 7.0914}, with corresponding eigenvectors 



1 



0.6814 0.4379 0.5865 



0.3883 -0.9048 0.1747 



Since £(£7, 0.04, 0.08)^ =v 1} by Corollary [2] and ([29} the eigenpair of M a> p is 



Ai = Ai + 



si 



2 + 1 = 3, = v\ . 



Moreover, £(£7,0.04,0.08)^ ^ v 2 , andbydMD, (A-I)' 



0.222 -0.0556 -0.3889 
-0.0556 0.1389 -0.0278 
-0.3889 -0.0278 0.8056 



[(A-I)Y = 
approximations 



0.2037 -0.0093 -0.3981 
-0.0093 0.0231 -0.0046 
-0.3981 -0.0046 0.8009 



A 2 = 1.5890, v\ 



. Then, (|251) and (1261) lead to the second-order 



0.7031 0.4409 0.5579 
in contrast to the first-order approximations of the eigenpair of M a g, which are equal to 



A 2 = 1.5933, vi 



0.6257 0.4242 0.6546 
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Moreover, £(£,0.04,0.08)^3 / v 3 , and by ([MI), (A-7I)~ 



-0.1511 0.0556 0.0244 
-0.0556 -0.0278 -0.0278 
0.0244 -0.0278 -0.1878 



[(A - 7/) 2 ]+ = 
proximations give 



0.0265 0.0093 -0.0067 
0.0093 0.0046 0.0046 
-0.0067 0.0046 0.0366 



Also, by (|2"5j) and (|2"6"j) . the second-order ap- 



A 3 = 7.0910, vi 



0.3878 -0.9049 0.1754 



but the first-order approximations are not satisfactory enough 



A 3 = 7.0867, vi 



0.3851 -0.9056 0.1779 
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